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Abstract. The neutron-proton mass difference in (isospin asymmetric) nuclear matter and finite nuclei is stud- 
ied in the framework of a medium-modified Skyrme model. The proposed effective Lagrangian incorporates both 
the medium influence of the surrounding nuclear environment on the single nucleon properties and an explicit 
isospin-breaking effect in the mesonic sector. Energy-dependent charged and neutral pion optical potentials in 
the s- and p-wave channels are included as well. The present approach predicts that the neutron-proton mass 
difference is mainly dictated by its strong part and that it markedly decreases in neutron matter. Furthermore, 
the possible interplay between the effective nucleon mass in finite nuclei and the Nolen-Schiffer anomaly is dis- 
cussed. In particular, we find that a correct description of the properties of mirror nuclei leads to a stringent 
restriction of possible modifications of the nucleon's effective mass in nuclei. 



1 Introduction 

This presentation is about a calculation of the neutron- 
proton mass difference in finite nuclei and in nuclear matter 
in the framework of a medium-modified Skyrme model. 

The nuclear-density dependence of isospin-breaking ef- 
fects belongs to one of the fundamental questions in nu- 
clear physics II1I2I3I . In particular, the effective neutron- 
proton mass difference in finite nuclei and nuclear mat- 
ter Am* av , although it is of astrophysical relevance, of im- 
portance for the description of mirror nuclei, the stabil- 
ity of drip-line nuclei etc., is not very well under con- 
trol. Rather, there exist very different theoretical predic- 
tions of this quantity for isospin-asymmetric nuclear mater 
II4I5I6I7I8I9I10I1 1I12I13I14I15I16I17 1 which predict both 
qualitatively and quantitatively different results. 

Most notably, the in-medium neutron-proton mass dif- 
ference is relevant for the description of mirror nuclei and 
the Nolen-Schiffer anomaly (NS A) H18I191 . Although there 
are many theoretical approaches devoted to the explanation 
of the NSA discrepancy 12012 1 122123 124125 126I27I28I29I30I3 1 1 
this phenomenon is still not fully understood. 

Skyrme-soliton models have the inherent advantage, 
compared with other hadronic models, that they are based 
on chiral input (chiral Lagrangians in the meson sector) 
and that they treat the structure, properties and interac- 
tions of the investigated nucleons on an equal footing 
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032I33I34I35I36I37I . Moreover, Skyrme-soliton models al- 
lows for an easy incorporation of the nuclear background 
via local medium-dependent coefficients [38 39 40]. They 
are therefore well suited to address the problem of the me- 
dium-dependence of the effective neutron-proton mass split- 
ting and the Nolen-Schiffer anomaly. In Refs. [42. 43 1 we 
investigated Am* nv in symmetric and asymmetric nuclear 
matter, respectively. Here we rather focus on the results 
of Ref. 1 44 1 which addresses the behavior of this quantity 
in finite nuclei and the Nolen-Schiffer anomaly. 

The presentation is organized as follows. In Sect. [2] we 
review how the Skyrme model should be generalized in or- 
der to include explicit isospin breaking and explicit chiral 
breaking on the same footing. An unconventional break- 
ing mechanism is introduced and the relevant quantization 
procedure is discussed. Section[3]is dedicated to the intro- 
duction of medium-modifications to the isospin-extended 
Skyrme model. We review the construction of the medium- 
modified Skyrme model based on the pion dispersion in a 
nuclear background. Section|4]focusses on the static results 
of the nucleon when the pertinent skyrmion is placed at the 
center of spherical symmetric nuclear core. The construc- 
tion of the electromagnetic part of the in-medium neutron- 
proton mass splitting is indicated. Values for the effective 
proton mass, the (total and electromagnetic contribution of 
the) in-medium neutron-proton mass splitting, the proton 
and neutron magnetic moments and the scalar and isovec- 
tor rms-radii are reported. In Sect. [5] we discuss the in- 
medium neutron-proton mass difference when the nucleon 
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is located at a distance R from the center. Corresponding 
results in infinite nuclear matter can be found in Sect. [6] In 
Sect. [7] tentative conclusions are given to the behavior of 
the effective shift dm* np . The Nolen-Schiffer anomaly and 
the prediction of the medium-modified generalized Skyrme 
model are discussed in Sect. [8] We end the presentation 
with final remarks in Sect. [9] 



with m n ± = m n± denoting the strong-interaction part of 
the mass of the charged pions. Note that the electromag- 
netic contribution to the latter mass is beyond the frame- 
work of the present model and will not be considered here. 



2.2 Quantization 



2 Isospin-modified Skyrme model 
2.1 Generalized symmetry breaking 



The standard Skyrme model H32I33I34I35I36I37I consists 
of the non-linear <x-model Lagrangian, which is of sec- 
ond order in the derivatives, and a fourth order soliton- 
stabilizing term: 

p2 1 
L = if Tr ( a " t/3 " C/t ) + ^2 Tr[C/+ ^ C/ ' U^vUf, (1) 



32e 2 



where 



U(x) = exp(2i Y TV(x)/F n ) = 1 + i^Ll* 



(2) 



is the usual chiral SU(2) matrix formulated in terms of 
the pion field n" . The pion-decay constant F n (its physical 
value would have been 186 MeV) and the stabilizing pa- 
rameter e are normally - and also here - adjusted in such 
a way, that, after rigid-rotator quantization in the (baryon 
number) B- \ soliton sector, the empirical isospin-averaged 
nucleon and delta masses are fitted. Note that the usual 
Skyrme model preserves isospin symmetry. This holds even 
if there is a chiral symmetry breaking term 



£sB = ^Tr(C7 + C/t_ 2 ) 



(3) 



present, which induces the pion-mass term in the (B=0) 
meson-sector of the model, see especially [33 35 3 61371 . 
In order to incorporate explicit chiral symmetry breaking 
and isospin-breaking on equal footing as it is the case for 
the free pion Lagrangian 

Xmes = d M n + d^n- - m 2 n± 7i + n- + \ (d^&n - m*,>rV) , 



^=(7Ti + in 2 ) , 



the original Skyrme model has to be modified. Following 
Rathske PTD we therefore add to the ordinary Skyrme La- 
grangian ([1) the following generalized symmetry breaking 
term 



srSB 



F 2 I 
-Tr[(t/ 



(U-l)M (I/T-l) 
- l)T 3 M 2 A (U f - l)r 3 ] 



(4) 



which explicitly breaks chiral symmetry and isospin sym- 
metry. Here the following definitions have been applied: 

Ai= tjml+mlc), M 2 A = \(m\ ± -m 2 ^ (5) 



Following Refs. [42 43 44 1, we quantized the isospin-modified 
Skyrme model with the help of one constrained {tpi (t)) and 
three unconstrained time-dependent collective coordinates 

<jpi(t),<p 2 (t),<p 3 (i)) 

U At/A 1 ' with A -> A(f) = exp (i r • <p/2J (6) 
where 



<p\ = (1)1, (pi - a>2, (p-i = W3 + tp 3 = U>3 + a c 



(7) 



are the pertinent angular velocities. Under the spherical 
hedgehog ansatz t/(r) = exp (ir • rF(r)), there exist also 
collective coordinates linked to static rotations of the skyr- 
mion in coordinate space, where Q\ , Q 2 , Q 3 are the relevant 
angular velocities. In terms of both classes of angular ve- 
locities the spatially integrated Skyrme Lagrangian (Q~|i in- 
cluding the generalized symmetry breaking term (|4]i turns 
into 



-C2 + -C4 + -C 



= -M NP - M A A 2 



+ i(A 2 + A 4 ) 



(o)[-Qi) 2 + (a) 2 -Q) 2 + ((Di-Qi+a c ) 2 



= -M NP + \A Xti (w; - Od 2 + A (w 3 - Q 3 ) a c , 



(8) 



where Mnp is the static soliton mass when it is non-perturbed 
(NP) by isospin-breaking. It is the hedgehog mass in case 
the above-mentioned hedgehog ansatz is inserted for the 
time-independent chiral matrix. A 2 and A4 are the moments 
of inertia of the non-linear <x model Lagrangian £ 2 and the 
fourth-order stabilizing term £.4, respectively. In the last 
line of Eq. (O, we employed the constraint 



MaA 2 = \Aa 2 c 



(9) 



whose role is to balance the isospin-breaking terms, such 
that the hedgehog ansatz of the usual Skyrme model can 
still be applied. In terms of the canonical conjugate mo- 
mentum operators T, = dL/dui (and 7,- = dL/dQi = -T,) 
corresponding to isospin (and spin), the pertinent Hamilto- 
nian becomes 

H = M NP + \Aa] + g - a c f 3 , 

while the baryon states are \T, T3, J -T, J 3 = -T3). When 
the Hamiltonian is sandwiched by these states, the energy 
of a nucleon of isospin component : 



: \ is determined as 



Wp/n = Ep/ D - M N p + \Aa 2 c + ^ + a c \ 



(10) 
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The latter equation allows to isolate the strong part of the 
neutron-proton mass difference as zlmL 8 = a c with 



a c = 



2M A 
A 



(1.291 +0.686) Me V. (11) 



For consistency, the empirical electromagnetic contribu- 
tion -0.686 MeV to the in-vacuum value 1 .291 MeV of the 
neutron-proton mass difference had to be subtracted. 



3 Medium-modifications 

The explicit isospin-breaking term is only one part con- 
tributing to the neutron-proton mass difference in a nu- 
cleus or nuclear matter. The other part comes from the 
isospin asymmetry of the nuclear background. In order to 
include this asymmetry into the nuclear-background pa- 
rameters of a medium-modified Skyrme model, used in 
Refs. [38 39 40|, we start out by introducing medium-modi- 
fications to the free pion Lagrangian ©. 

Note that throughout this presentation we will denote 
medium-modified Langrange terms, energies, masses, mo- 
ments of inertia, form factors, etc. with an explicit asterix. 



3.1 Pion dispersion in a nuclear background 

The pertinent medium-modified Lagrangian of the pion field 
can be formulated in terms of the self-energies (i.e. energy- 
dependent optical potentials) of the charged pions as fol- 
lows: 

A=±,0 

= -Cmes - ij^fl ^ 2" 71 a + XS abi^a ~ ^ Kfrj. (12) 

In the local density approximation the i-wave self-energies 
become 

nf(aj,r) = -tej^r);*!^^ (13) 

with the total and isosvector density 

p(r) = p„(r) + p p (r) , dp(r) = p„(r) - p p (r) (14) 

and 77 = 1 + m„/mN in terms of the isospin-averaged pion 
and nucleon masses. The effective isoscalar and isovector 
pion-nucleon scattering lengths can be expressed as 

bf(aj) « bo(oj) - ^[b 2 (oj) + 2b\{uj\ , (15) 

bo(aj) * -b (l - m- 2 a> 2 Y^ri) , (16) 

bi(co) ~ hlm^co + 0.143m; V )/(4nrf) (17) 

in terms of the total Fermi momentum Icf and the parame- 
ters bo - -1.206 rrC K l ,b\ = —1.115 m~ l . The corresponding 



p-wave self-energies read in the local density approxima- 
tion 

fr+, , „ 47rc ± (w,r) „ 
p 1 + Ang'c^o, r) 



(V 2 c>,r)) (18) 



with g' = 0.47 and 

c*(w,r) = (coMpW + dMw)//, (19) 

in terms of the isoscalar and isovector pion-nucleon scat- 
tering volumes cq = 0.21m n 3 and c\ = 0.165 m~ 3 , respec- 
tively; see Refs. 0431441 for more details on the parameters 
and references. 

3.2 Medium-modified generalized Skyrme model 

Marrying the generalized isospin-modified Skyrme model 
of Rathske BP with the in-medium modified Skyrme mo- 
del of Refs. [38 39 40| (which is characterized by den- 
sity dependent coefficients of the second order symmetry- 
conserving and breaking Lagrangian terms) and extend- 
ing it to asymmetric nuclear background as in Eq. (TT2l . 
one arrives at the following isospin- and medium-modified 
Skyrme Lagrangian M42I43I441 : 



X." = t\ + L\ + -£* S B + ^-Cmes + 4-C en 



(20) 



^=M l+ %) TT{doUdoUf) 

-(l-aQTrfVtf.VI/t)}, 



£4 = 



1 



32e 2 



r n V" 1 

16 2-i 

a=l 



Tr(T fl f/)Tr(T fl f/ T ), 



a,b=\ 



2m„ 



Tv(T a U)Tv(T b d tf) . 



Here £, mes is the explicitly isospin-breaking Lagrangian, 
while ^f-Cenv summarizes the enviroment-induced isospin- 
breaking contribution. The second-order Lagrangians £f 2 
and XfSB* are isospin-symmetric, but explicitly medium- 
dependent, whereas X4 is neither isospin-breaking nor me- 
dium-dependent. We preserve the original term for stabi- 
lizing skyrmions and do not modify it by, e.g. higher-order 
derivative terms or vector-meson contributions, since it is 
on the one hand generic for our purposes and on the other 
hand still simple to handle. 

The medium functionals appearing in £1 and -C^-sb are 
constructed by the rule 77(w) — > n(ido) from the input of 
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Sect. 13. II and take the form: 



3k f 



4 = 



8tt 2 (1 + m n jm^) 
3k F 

An 2 (l + to^/toai) 
2nc + 2nc~ 



*o P. 



1 + 4ng'c + 1 + \ng'c- 



with 



± _ cop + ci^p 
1 + m K jm^ ' 

4^ = £ 1( 5p, Ax P 



2mn„ 



(21) 
(22) 
(23) 

(24) 

Cl (v 2 <5p). (25) 



Note the explicit occurrence of gradient density terms in 
the last expression which especially operate at the surface 
of finite nuclei and vanish in homogenous nuclear matter; 
compare Ref. flU) with l43l . 



3.3 Quantization of the in-medium Skyrme model 



By essentially copying the steps of Sect. 12.21 one can find 
the following expression for the spatially integrated medium- 
dependent Lagrangian d20l >, see Refs. |43 44 1 for more de- 
tails: 



d 3 r£* 



-Ml 



r. 0/7 + 

2 A 4- 1 2 A* 

■j-'imes ' n /1 lolo,12 



(a)iQi + cj2Qz)Al )Ql2 
(a*i - Qt, + a*) 2 



-A 



wfl,33 



Q 2 +Ql 

9 n QQ,\2 



+ (a>3 - Q?, + fl*)A* n 



Here w, and £?, are the angular velocities in isotopic and co- 
ordinate space, respectively, while A* w .., Al QJj and A* QQij 
are the pertinent medium-dependent moments-of-inertia. 
Finally, A* nv is the medium-dependent coefficient of the an- 
gular velocity sum (lot, —Q$ +a*) generated by the environ- 
ment-induced isospin breaking term A£* em . 

Constructing the pertinent Hamiltonian by standard means 
and sandwiching it between the usual isospin-spin states, 
one can determined the energy of a nucleon (or delta) state 
classified by the total isospin T and third component Tt, as 



E = M* NP + VM 2 A mes 



A" 1 

1 *env 

A* + A* _ 2/\* 
/l QO,12 /l bib),12 ^ ll iaQ,Yl 



A* - A* 2 ) 



2(a: 



ojoj, 12 QQ, 12 



1 3 



2 A* 



a,, + 



Aenv 
^wfi,33 J 



(t(t + 1) - r 2 ) 



T 3 . 



(26) 



Consequently, the strong part of the neutron-proton mass 
difference in the interior of a nucleus can be isolated from 
the last term of Eq. d26l i and takes the form 



A *(strons) * 

Aninp 6 = a c + 



A* 



(27) 



with a] 



= ^M 2 A A n 



;s /^wQ,33- 



4 Skyrmion at the center of a nucleus 



The knowledge of Eqs. d26l i and d27l i allows to construct 
the total mass of a proton (neutron) and the strong contri- 
bution to neutron-proton mass difference when the nucleon 
is located at the core of a nucleus. 

In order to determined the total in-medium neutron- 
proton mass splitting, the electric (E) and magnetic (M) 
form factors of the nucleons must be utilized in addition: 

Ge(Q 2 ) = £ Jd 3 r***Ar), 
Gm(Q 2 ) = ^ J dW^rxKr)], 

where q is the transferred momentum. Furthermore, j° and 
j correspond to the time and space components of the prop- 
erly normalized sum of the density-dependent baryonic cur- 
rent B* and the third component of the isovector current V* 
of the Skyrme model; for more details see Refs. [43 44]|- 

The electromagnetic part of the neutron-proton mass 
difference is then calculated in terms of the scalar (S) and 
isovector (V) electric (E) and magnetic (M) form factors 
according to the formula 



Am 



t(EM) 
np 



CO 



^2-^(q 2 )G£(q 2 ) 



(28) 



of Ref. B31 . where a is the electric fine structure constant. 

Explicit expressions of the pertinent scalar and isovec- 
tor charge densities (Pg v ) and magnetic densities (p^ V ) 
can be found in Ref. |44). The rms radii (r 2 )* 1 /, 2 and (r 2 )* 1 ^ 2 
are calculated from the charge densities in the standard 
way, whereas the in-medium magnetic moments of the pro- 
ton and neutron follow from calculating the sum and differ- 
ence of the integrated scalar and isovector magnetic den- 
sities, respectively. Explicit expressions for the in-medium 
skyrmion mass M^ p and the moments-of-inertia appearing 
in Eqs. (|26| | and d27l i can be found in Refs. [44] as well. 

In summary, the (in-medium) mass of the proton, the 
total and electromagnetic contribution of the (in-medium) 
neutron-proton mass splitting, the (in-medium) magnetic 
moments of the proton and neutron and the isocalar and 
isovector (in-medium) rms radii are listed in Table Q] for 
the following situation: either a skyrmion in free space or 
located at the center of a finite nucleus core (e.g. 14 N, le O, 
38 N, 4() Ca) is quantized as a proton (or neutron). In the case 
that it is quantized as a proton, the resulting total nuclei 
will be I5 , 17 F, 39 Ca, 41 Sc, respectively; in the case of 
a neutron, the resulting nuclei are 15 N, n O, 39 N, 41 Ca in- 
stead. 

In accordance with previous calculations for infinite 
nuclear matter in the isospin-symmetric case [42|, the sec- 
ond and third columns of Table [TJ show that the total and 
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Table 1. Static properties of a nucleon when the pertinent skyrmion is either in free space of added to the center of a finite nucleus 
core, namely 14 N, 16 0, 38 K, and 40 Ca, respectively, such that a proton or the total nucleus 15 0, 17 F, 39 Ca, and 41 Sc results. Here m* is the 
in-medium proton mass, /fm* p is the in-medium neutron-proton mass difference and Am^ p EM} is its electromagnetic part. Furthermore, yU* 
and yu* are the in-medium proton and neutron magnetic moments in units of the free-space Bohr magnetons (n.m.). Finally, (r 2 )^ 2 and 
(r 2 )*'y 2 are the in-medium isoscalar (S) and isovector (V) charge radii of the nucleon inside a nucleus. For more details see Ref. 1441 . 



Core->Elem. nf p Am' np Am™ ^ £ (r%f (r 2 )^ 

[MeV] [MeV] [MeV] [n.m.] [n.m.] [fm] [fm] 



Vac— » p 


938.268 


1.291 


-0.686 


1.963 


-1.236 


0.481 


0.739 


M N ^ 15Q 


593.285 


1.668 


-0.526 


2.355 


-1.276 


0.656 


0.850 


16 0-> 17 F 


585.487 


1.697 


-0.517 


2.393 


-1.297 


0.667 


0.863 


38 K -> 39 Ca 


558.088 


1.804 


-0.480 


2.584 


-1.422 


0.722 


0.942 


40 Ca-> 41 Sc 


557.621 


1.804 


-0.478 


2.569 


-1.428 


0.724 


0.947 



electromagnetic part, respectively, of the neutron-proton 
mass difference are slightly increased in finite nuclei rel- 
ative to the vacuum case. Furthermore, the present model 
predicts (see the first column of Table [TJ that the effective 
nucleon mass is strongly reduced at the center of the nu- 
cleus. We will return to this point later on when we discuss 
the Nolen-Schiffer anomaly in Sect. [8] 

Finally, the increase in the tabulated in-medium val- 
ues of the magnetic moments and rms charge radii is com- 
patible with the hypothesis of the nucleon-swelling in a 
nuclear background which was already confirmed by the 
isospin-symmetric calculations of Refs. [38 39 40] for the 
medium-modified Skyrme model; see also the static in-me- 
dium results of Refs. 1421431 . 



5 Skyrmion off-centered in a nucleus 



Fig. 1. A sketch of a skyrmion located inside a finite nucleus with 
R = |R| the separation between the geometrical center (O') of the 
skyrmion and the center (O) of the nucleus. The vectors (angles) 
r (6) and r' (9') refer to the body-fixed coordinates of the nucleus 
and skyrmion, respectively. Since the nucleus is spherically sym- 
metric, both coordinate systems can be orientated in such a way 
that their z-axes coincide. 





V.9' \ 




1 o Rl 


O' J 


z 



Whereas the results of the last section hold for the case 
that the skyrmion is at the center of a symmetrical core 
of a nucleus, the results of this section apply for the more 
general situation that the skyrmion is located at a distance 
R from the center of a finite nucleus, as specified in Fig.[T] 
Note that in this case the spherical symmetric hedgehog 
ansatz cannot be used anymore, since the background - as 
viewed from the center of the skyrmion - is not spherically 
symmetric. Thus there exists further deformations in the 
isotopic and in the coordinate space [38 39]. In this case a 
variational ansatz can still be used which allows for more 
freedom in the radial and especially in the angular coordi- 
nates than the hedgehog form would do, but is still com- 
patible with the quantization procedure. See Ref. fl4ll for 
more details on this variational parameterization. In this 
reference, it was also checked that the results of the vari- 
ational computation merge with the ones of the hedgehog 
computation when the skyrmion is moved back to the cen- 
ter of the nucleus. 



between the geometrical center of the pertinent skyrmion 
and the center of the core-nucleus, i.e. 14 N, le O, 38 K, and 
40 Ca. When the skyrmion is quantized as a proton, the re- 
sulting nuclei are ls O, 17 F, 39 Ca, and 41 Sc, respectively. 
Note that the in-medium mass of the proton m* starts with 
the value listed in Table[T]when its is near the center of the 
nucleus - in agreement with the statements at the end of 
the last section. With increasing distance from the center 
of the nucleus the value of the effective mass monotoni- 
cally increases, until it smoothly approaches the free space 
value m p at the border of the nucleus. 



5.1 Effective proton mass m* inside a nucleus 

First, we will show in Fig. [2] the results of Ref. [44] for 
the effective proton mass to* as function of the distances 



We will come back to the effective nucleon mass when 
we discuss the Nolen-Schiffer anomaly in Sect. [8] In the 
following we will concentrate on the effective in-medium 
neutron-proton mass difference. 
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Fig. 2. The dependence of the effective mass of the proton m* 
on the distance R between the center of the skyrmion and the 
center of the nucleus. The solid curve represents the case that a 
skyrmion - quantized as a proton - is added to a 14 N nucleus, 
giving 15 O in total. The dot-dashed curve refers to the case of a 
16 nucleus — > 17 F in total. The dashed curve stands for the case 
of 38 K — > 3, Ca and the dotted curve represents the case of 40 Ca — » 
41 Sc. The horizontal line marks the free space value of the proton 
mass. Figure from Ref. |44|. 

1000 1 i 




8 



R[fm] 



5.2 Effective neutron-proton mass difference 

In Fig. [3] the behavior of the strong part of the in-medium 
neutron-proton mass difference can be found for the same 
notations and input as in Fig. [2] 



clear background is high, the value of Am„ p is at an 
extremum. 

In Fig. 2] the electromagnetic part of the in-medium 
neutron-proton mass difference is shown for the same nu- 
clei as before. It can be seen that the variations in the elec- 
tromagnetic part of the effective neutron-proton mass dif- 
ferences are small as compared with their strong counter 
parts. Furthermore, z)m*p EM) nearly monotonically decreases 
with increasing distance from the center. Note that the val- 
ues at R = agree of course with the values listed in Ta- 
ble [Q 

Fig. 4. The dependence of Am^ M> on the distance R between 
the center of the skyrmion and the center of the nucleus. The 
notations and input are the same as in Fig. [2] with the excep- 
tion that the horizontal line marks the free space value Am'™ ) = 
-0.69 MeV. Figure from Ref. |44) . 



-0.4 r 




4 

R ffml 



Fig. 3. The dependence of Am^ Uons) on the distance R between 
the center of the skyrmion and the center of the nucleus. The 
notations and input are the same as in Fig. [2] with the exception 
that the horizontal line marks the free space value Am^° ng) = 
2.0 MeV. Figure from Ref. (44). 



Finally, for completeness, in Fig. [5] the dependence of 
the total effective neutron-proton mass difference Am* on 
the distance R is shown for the same nuclei as specified 
above. Again, the values at R = agree with the values 
listed in Table Q] 



2.4 r 



> 




2 4 6 

R [fm] 



Note that the strong part of the in-medium mass dif- 
ference has a non-monotonic behavior. This follows from 
the fact that the density is a local quantity and that addi- 
tional isospin-breaking contributions arise due to the den- 
sity gradients resulting from the p-wave pion-nucleus scat- 
tering, see Eq. (t25t . Especially, it can be observed that in 
the surface region of each nucleus, where the density gra- 
dients are large and the local isospin asymmetry in the nu- 



Fig. 5. The dependence of the total in-medium neutron-proton 
mass difference Am* nv on the distance R between the center of the 
skyrmion and the center of the nucleus. The notations and input 
are the same as in Fig. [2] with the exception that the horizontal 
line marks the free space value Am np = 1.3 MeV. Figure from 
Ref. (441 . 
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6 Neutron-proton mass difference in 
nuclear matter 

The densities of the finite nuclei discussed here and in 
Ref. [44 1 correspond closely to the isosymmetric nuclear 
matter case studied in Refs. [4243], where a moderate in- 
crease of the effective neutron-proton difference, compati- 
ble with Fig.[3]or Fig. [5] can be observed. In neutron mat- 
ter, however, as Ref. Il43l shows, there is a pronounced de- 
crease of the in-medium neutron-proton mass difference 
with increasing density. 

Therefore, instead of discussing the density-variations 
of all the static quantities in an infinite nuclear matter back- 
ground as in Ref. Il43l . we report here only about the ef- 
fective neutron-proton mass differences in nuclear matter 
studied in Ref. (431 . 

In Fig. [6] the strong part of the in-medium neutron- 
proton mass splitting in nuclear matter, Am* , is shown for 
isospin-symmetric nuclear matter (solid curve), neutron- 
rich matter (dashed curve), pure neutron matter (dotted cur- 
ve), and proton-rich matter (dot-dashed curve). Especially, 
when the isospin symmetry of nuclear matter is broken, 
z)m*p tr strongly varies (see the dashed and dot-dashed 
curves in Fig. |6]). In pure neutron matter the change be- 
comes very drastic (see the dotted curve in Fig. |6), and 
/lm np strong) decreases very rapidly with increasing density. 



Fig. 6. Density dependence of the strong part Arri!^ of the 
neutron-proton mass difference in nuclear matter. The abscissa 
represents the density p normalized to the saturation density of 
ordinary nuclear matter po = 0.5ml, while the ordinate shows the 
mass difference in units of MeV. The result in isospin-symmetric 
matter is plotted as a solid curve, the result of neutron-rich matter 
with Sp/p = 0.2 as dashed curve, the dotted curve represents pure 
neutron matter (Sp/p = 1) and the dot-dashed curve proton-rich 
matter with Sp/p = -0.2. Figure from Ref. 1 43 1. 
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In contrast to the strong part, the electromagnetic part 
of the neutron-proton mass difference in nuclear matter 
varies only by a very small amount when the isospin-asym- 
metry parameter 6p/p is increased (see Fig . O . 

Again, for completeness, the total neutron-proton mass 
difference in nuclear matter is presented in Fig. [8] The dif- 
ference to the purely strong case, shown in Fig. [6] is hardly 
visible. 



Fig. 7. Density dependence of the electromagnetic part Am*J- p ' 
of the neutron-proton mass difference in nuclear matter. The axes 
and curves are defined as in Fig. [6] Figure from Ref. 1431 . 
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Fig. 8. Density dependence of the total neutron-proton mass dif- 
ference Am* in nuclear matter. The axes and curves are defined 
as in Fig. [6] Figure from Ref. [43 1. 
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7 Tentative Conclusions 

In summary, we have studied the effective neutron-proton 
mass difference Am* np in finite nuclei in the framework of 
an isospin- and medium-modified Skyrme model. 

The in-medium mass of the proton starts near the cen- 
ter of the nucleus with a strongly reduced value and in- 
creases monotonically with increasing distance from the 
center. The strong part of the effective neutron-proton mass 
difference has a non-monotonic behavior because of the 
additional isospin-breaking contributions due to the gradi- 
ents arising from the p-wave pion-nucleus scattering. There 
is an extremum of A™ 10118 at the surface of the nucleus be- 
cause there the gradients are large and the isospin asym- 
metry in the nuclear background is high. The electromag- 
netic part of the effective neutron-proton mass difference 
is negative and increases monotonically in magnitude with 
increasing distance from the center. In magnitude its value 
is small as compared with its strong counterpart. The other 
static quantities (in-medium magnetic moments, rms charge 
radii etc.) behave according to the nucleon-swelling hy- 
pothesis for a nucleon embedded in the medium. 

While in isosymmetric backgrounds the effective neu- 
tron-proton mass difference moderately increases, in asym- 
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metric nuclear matter the strong part of this quantity can 
vary markedly (upwards in proton-rich matter and down- 
wards in neutron-rich matter), whereas the electromagnetic 
contribution is always small and subleading. Especially in 
neutron matter, there is a strong decrease of the (strong and 
also the total) neutron-proton mass difference with increas- 
ing density. 

8 Nolen-Schiffer anomaly 

A long standing problem in nuclear physics is the Nolen- 
Schiffer anomaly observed in mirror nuclei II 1 81 1 91 . Here 
we will show how the Nolen-Schiffer anomaly can be treated 
in the framework of an isospin- and medium modified gen- 
eralized Skyrme model. 

The mass difference between mirror nuclei 

AM = z+l M N - A M N+1 (29) 

which differ by one unit in their charges, AZ = 1, can very 
precisely be measured and is usually split into two terms 

AM = AE EM - Am np - Liml p - Am^J . (30) 

The first contribution is the Coulomb energy difference 
AEem computed relatively to the free neutron-proton mass 
difference, whereas the second is the in-medium neutron- 
proton mass difference subtracted from the free one and 
therefore the net in-medium change of this quantity. The 
Coulomb energy difference includes various corrections, 
e.g., by exchange terms, the center-of-mass motion, finite 
size effects of the proton and neutron distributions, mag- 
netic interactions, vacuum effects, short-range two-body 
correlations etc. It can be calculated with great accuracy 
(within 1 % error) |[T9l . If Am* np is assumed to be constant 
and equal to the vacuum value, then Eq. (l30l cannot be 
satisfied. This phenomenon is called the Nolen-Schiffer- 
anomaly (NSA). Quantitatively, the NSA ranges - through- 
out the periodic table - from a few hundred keV for the 
lightest nuclei up to several MeV for the heaviest ones. 
A possible resolution is the assumption that the effective 
neutron-proton mass difference would decrease with in- 
creasing mass number A, such that 

Ansa = Am np - Am* np . (31) 

Within the present approach Am* np has a local R depen- 
dence according to the location of the nucleons inside the 
nuclei as shown in Fig. [5] In order to compare these results 
with the experimental data one therefore has to average the 
value of Am* np with respect to the separation R. Since the 
nucleons inducing the Nolen-Schiffer anomaly are valence 
nucleons, these must be located in the peripheral region of 
each of the mirror nuclei, if the latter differ by one particle 
or hole from a (magic) closed-shell nucleus. The averaged 
effective masses and splittings can therefore be expressed 
as follows: 

m* a = f ml(R)\if, n (Rfd 3 R, 



m; = j m ;(R) |<A P (fl)| 2 d 3 fl , 

Am* np = ln a -ln p , (32) 

where \ifr a (R)\ 2 and \tfr p (R)\ 2 are the density distributions of 
the in-medium neutron and proton, respectively. In terms 
of the difference of the density distributions 

Aifi p (R) s I^Wl 2 - \MR)\ 2 . 

Eq. ( 1321 can be rewritten as 

Am np » J {Aifr 2 np (R) n$(R) + Am^R) \if, p (R)\ 2 } d 3 R 

= Am^+Am^\ (33) 
where the subleading contribution of the cross term 

f ^l^K^R 

is neglected. Thus here the Nolen-Schiffer anomaly simply 
reads 

zInsa = Am ap - \Am^ + Am*^ . (34) 

It is listed in Table |2j left panel (labeled 'a=0') of the 
third column (labeled 'Present approach') for chosen pairs 
of mirror nuclei, such that the overall mass number A in- 
creases. 

Whereas qualitatively the calculated NSA values have 
the correct A behavior, quantitatively the results are more 
than one order of magnitude too big (compare with the em- 
pirical values listed in the last column of Table[2]). This can 
be traced back to the pronounced negative shift of AnQ. 
(see the first entry of the third column of Table |2j. This 
shift is mainly there for three reasons: 

(i) the rather large renormalization of the effective nucleon 
mass, 

(ii) the pronounced R dependence of m* inside the nucleus 
(see Fig.|2]i, and 

(iii) the relative swelling of the proton distributions due to 
the Coulomb factor, i.e. Ai[> 2 w + 0. 

For example, the averaged in-medium mass of the valence 
proton in n O is reduced to m* p = 812.35 MeV. This drop 
of about 125 MeV is very large in comparison with the em- 
pirical value of the binding energy per nucleon in nuclear 
matter. For heavier nuclei, where the density in the interior 
approximates the normal nuclear matter density, the drop 
of the averaged effective mass is even larger, e.g. m p -m* p ~ 
(150 - 200) MeV in the 40 Ca region (see the second col- 
umn of Table |2]i down to ~ 300 MeV in the 208 Pb region. 

If solely the contribution Am^ 2) (due to the explicit R 
dependence of the neutron-proton mass difference) were 
considered, then the NSA in the present approach would 
even have a negative sign: Am np - Am*^ < 0. 

Instead of driving the input parameters of Sect. 13 . 1 1 to 
unphysical values in order to match the NSA discrepancy, 
we suggested in Ref. [44] to invert the problem and to 
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Table 2. The averaged mass m* of the valence proton in a given nucleus, the contributions to the effective neutron-proton mass difference 
(see Eq. d33t ) and the Nolen-Schiffer discrepancy zl NSA calculated in the present approach by Eq. i34\ or Eq. J35t and the corresponding 
"empirical" results of Ref. 1191 . All quantities are in units of MeV. 



Present approach 



Nuclei 




Q-rcn = 


«,,„ = 0.95 


— empirical 
^NSA 

Ref. fT9) 


«,cn=0 a rcn =0.95 


^NSA 


fiSg A ™t ^nsa 


ls O- 15 N 
n F .n 

39 Ca- 3< >K 
41 Sc- 41 Ca 


767.45 928.30 
812.35 930.54 
724.78 926.16 
771.71 928.51 


-4.27 1.56 4.02 
-5.53 1.52 5.33 
-8.11 1.67 7.75 
-9.74 1.62 9.44 


-0.21 1.33 0.20 
-0.28 1.32 0.27 
-0.41 1.33 0.37 
-0.49 1.33 0.47 


0.16 ±0.04 
0.31 ±0.04 
0.22 ± 0.08 
0.59 ± 0.08 



estimate the effective nucleon mass inside finite nuclei ac- 
cording to the results in the isospin-breaking sector. To per- 
form this task an artificially added renormalization param- 
eter a Kn in the expression 

m* np (R, a ren ) = m* p (R) + {m^ p - m* , p (R))a Ka (35) 

of the effective nucleon mass is fine-tuned in such a way, 
that the Nolen-Schiffer anomaly is satisfied. The results 
are presented in Table [2] fourth column labeled 'a ren = 
0.95'. It can be seen that a successful description of the 
correct order of the NSA implies a rather small drop of the 
mass of the valence nucleons: m ap - m* lp {a len = 0.95) ~ 
10 MeV which is close to the empirical binding energy per 
nucleon in nuclear matter. In this case, the contribution to 
the Nolen-Schiffer anomaly from the term Am"^ can be 

neglected: Am np - Am^ p 2 \a iea = 0.95) ~ -0.02 MeV. 

9 Final remarks 

Let us discuss the relevance of our results for the Nolen- 
Schiffer anomaly. Qualitatively, our approach predicts the 
sign and the relative mass-number increase of this anomaly. 
But quantitatively it is far from satisfactory: the results 
are more than one order of magnitude too large. Clearly, 
the part of our calculation relevant to the Nolen-Schiffer 
anomaly depends on the proton and neutron distributions 
of the mirror nuclei and is very sensitive to the behavior 
of the wave functions of the valence nucleons in the pe- 
ripheral region of the nucleus. We have pointed out the 
possibility that the Nolen-Schiffer anomaly may rather fol- 
low from the behavior of the effective nucleon mass in 
finite nuclei than from the effective neutron-proton mass 
difference: our calculations imply that the Nolen-Schiffer 
anomaly could not and, maybe, should not be saturated by 
z)m*p 2) (the averaged contribution due to the explicit den- 
sity and radial dependence of the neutron-proton mass dif- 
ference). Rather more important is Am*J p \ the contribu- 
tion due to the difference in the squared wave functions 
of valence proton and neutron weighted by the local (den- 
sity and density-gradient induced) variation of the effec- 
tive mass of the nucleon. In fact, when we restrict the in- 
medium reduction of the (averaged effective) proton mass 
to about 1 % of the free proton mass - a value which is 



compatible with the empirical binding energy per nucleon 
in nuclear matter - we obtain a rather precise description of 
the NSA. Here, we should remark that the gradient terms 
which are present in our model do not noticeably affect the 
scaling behavior of m*. They are important for the surface 
behavior of Am* np , though. 

The calculated Coulomb energy differences AEem d30l > 
for mirror nuclei of Ref. [19] incorporate a contribution 
due to the different wave functions of valence nucleons that 
is known as Thomas-Ehrman effect 1 46 47 1 . Note, how- 
ever, that in Ref. |[T9l a constant, /^-independent value of 
the nucleon mass - namely the free mass in vacuum - was 
used, whereas here the effect is based on the R-dependence 
of the effective nucleon mass. 

In summary, the possibility exists that the anomaly of 
the mirror nuclei can be saturated by invoking a dynamical 
(local) mass of the nucleon that needs to be only slightly 
reduced in comparison to its free counter part. In this con- 
text, it should be pointed out that the anomaly could also 
be removed by the introduction of an additional term into 
the Coulomb part of the energy density of nuclear systems 
in the local-density-functional approach [48 49 50 50- In 
fact, this term is chosen such that it is proportional to the 
isoscalar rather than to the isovector density. Furthermore, 
it should be stressed that the isoscalar contribution is sur- 
face dominated. It is even argued in Refs. [52 53 1 that the 
effective-to-free-nucleon-mass ratio m* N /niN is unity to within 
a few percent. Apparently, different model calculations can 
lead to similar conclusions about the origin of the Nolen- 
Schiffer anomaly. 

Returning to the in-medium modified Skyrme model, 
we would like to point out that the results might still be im- 
proved if the calculations can be made more self-consistent, 
e.g. , by the incorporation of feedback mechanisms between 
the modified skyrmion and the local nuclear background. 
In addition, the inclusion of further degrees of freedom 
might be a useful possibility as their introduction can any- 
how be motivated by more detailed considerations about 
the nucleon structure and the nucleon-nucleon interaction. 
Also the non-local character of the effective nucleon mass 
may be of importance. 



EPJ Web of Conferences 



Acknowledgement 

We would like to thank Horst Lenske for useful discus- 
sions and a computer code calculating nuclear densities of 
finite nuclei. We are also grateful to Frank Griimmer for 
providing us with calculated shell-model wave functions. 
The work of U.T.Y. was supported by the Alexander von 
Humboldt Foundation. The work of A.M.R. was supported 
by the second phase of the Brain Korea 21 Project in 2007 
and by the German Academic Exchange Service DAAD. 
Partial financial support from the EU Integrated Infras- 
tructure Initiative Hadron Physics Project (contract num- 
ber RII3-CT-2004-506078), by the DFG (TR 16, "Subnu- 
clear Structure of Matter") and by BMBF (research grant 
06BN411) is gratefully acknowledged. This work is par- 
tially supported by the Helmholtz Association through funds 
provided to the virtual institute "Spin and strong QCD" 
(VH-VI-231). 



References 

1. B.A. Li, CM Ko and W. Bauer, Int. J. Mod. Phys. E 
7 (1998) 147 CarXiv:nucl-th/97070141 

2. V. Baran, M. Colonna, V. Greco and M. Di Toro, Phys. 
Rept. 410 (2005) 335 |arXiv:nucl-th/0412060| | 

3. A.W. Steiner, M. Prakash, J. M. Lattimer and P.J. Ellis, 
Phys. Rept. 411 (2005) 325 [ arXi v:nucl-th704"T 0066 1 

4. M. Lopez-Quelle, S. Marcos, R. Niembro, A. Bouyssy 
and V.G Nguyen, Nucl. Phys. A 483 (1988) 479 

5. I. Bombaci and U. Lombardo, Phys. Rev. C 44 (1991) 
1892 

6. E. Chabanat, J. Meyer, P. Bonche, R. Schaeffer and 
P. Haensel, Nucl. Phys. A 627 (1997) 710 

7. S. Kubis and M. Kutschera, Phys. Lett. B 399 (1997) 
191 larXiv:astro-ph/97030491 

8. W. Zuo, I. Bombaci and U. Lombardo, Phys. Rev. C 
60 (1999) 024605 [ arXiv:nucl-th/ 102035] 

9. K. Tsushima, K. Saito and A. W. Thomas, Phys. Lett. 
B 465 (1999) 36 |arXiv:nucl-th/990710il 

10. V. Greco, M. Colonna, M. Di Toro, G. Fabbri and 
F. Matera, Phys. Rev. C 64 (2001) 045203 (2001) 

11. F. Hofmann, CM. Keil and H. Lenske, Phys. Rev. C 
64 (2001) 034314 (2001) |arXiv:nucl-tn70 007050| 

12. B. Liu, V. Greco, V. Baran, M. Colonna and M. Di 
Toro, Phys. Rev. C 65 (2002) 045201 

13. W. Zuo, L.G Cao, B.A. Li, U. Lombardo and 
C.W Shen, Phys. R ev. C 72 (2005) 014005 
|arXiv:nucl-th/0506003l 

14. E.N.E. van Dalen, C. Fuchs and A. Faessler, Phys. 
Rev. Lett. 95 (2005) 022302 ]arXiv:nucl-th/05020641 

15. T. Lesinski, K. Bennaceur, T. Duguet and 
J. Meyer, Phys. R ev. C 74 (2006) 044315 
|arXiv:nucl-th/0607065l 

16. E.N.E. van Dalen, C. Fuchs and A. Faessler, Eur. Phys. 
J. A 31 (2007) 29 |arXiv:nucl-tn706T2 066| 

17. L .W Chen, CM . Ko, B.A. Li and G.C. Yong, 
larXiv:0704.2340l [nucl-th] 

18. J.A. Nolen and J.P Schiffer, Ann. Rev. Nucl. Part. Sci. 
19(1969)471 



19. S. Shlomo, Rep. Prog. Phys. 41 (1978) 957 

20. S. Shlomo, Physica Scr. 26 (1982) 280 

21. E.M. Henley and G. Krein, Phys. Rev. Lett. 62 (1989) 
2586 

22. T. Hatsuda, H. Hogaasen and M. Prakash, Phys. Rev. 
Lett. 66 (1991) 2851 [Erratum-ibid. 69 (1992) 1290] 

23. A.G Williams and A.W. Thomas, Phys. Rev. C 33 
(1986) 1070 

24. T.D. Cohen, R.J. Furnstahl and M.K. Banerjee, Phys. 
Rev. C 43 (1991) 357 

25. M.H. Shahnas, Phys. Rev. C 50 (1994) 2346 

26. C.J. Horowitz, J. Piekarewicz, Phys. Rev. C 63 (2000) 
011303R 

27. T. Suzuki, H. Sagawa and A. Arima, Nucl. Phys. A536 
(1992) 141 

28. T. Schafer, V. Koch and G.E. Brown, Nucl. Phys. A 
562(1993)644 

29. E.G. Drukarev and M.G. Ryskin, Nucl. Phys. A 572 
(1994)560 

30. C. Adami and G.E. Brown, Z. Phys. A 340 (1991) 93 

31. B.K. Agrawal, T. SO, S.K. Samaddar, J.N. De and 
S. Shlomo, Phys. Rev. C 64 (2001) 024305 

32. T.H.R. Skyrme, Proc. Roy. Soc. Lond. A 260 (1961) 
127 

33. T.H.R. Skyrme, Nucl. Phys. 31 (1962) 556 

34. GS. Adkins, C.R. Nappi and E. Witten, Nucl. Phys. B 
228(1983)552 

35. G.S. Adkins and C.R. Nappi, Nucl. Phys. B 233 
(1984) 109 

36. I. Zahed and G.E. Brown, Phys. Rep. 142 (1986) 1 

37. U.-G. MeiBner and I. Zahed, Adv. Nucl. Phys. 17 
(1986) 143 

38. U.T Yakhshiev, M.M. Musakhanov, A.M. Rakhimov, 
U.-G. MeiBner, and A. Wirzba, Nucl. Phys. A700 
(2002) 403 

39. U.T. Yakhshiev, U.-G. MeiBner and A. Wirzba, Eur. 
Phys. J. A16 (2003) 569 

40. U.T. Yakhshiev, U.-G. MeiBner, A. Wirzba, 
A.M. Rakhimov, and M.M. Musakhanov, Phys. 
Rev. C71 (2005) 034007 

41. E. Rathske, Z. Phys. A 331 (1988) 499 

42. U.-G. MeiBner, A.M. Rakhimov, A. Wirzba, and 
U.T. Yakhshiev, Eur. Phys. J. A31 (2007) 357 

43. U.-G. MeiBner, A.M. Rakhimov, A. Wirzba, and 
U.T. Yakhshiev, Eur. Phys. J. A32 (2007) 299 

44. U.-G. MeiBner, A.M. Rakhimov, A. Wirzba, and 
U.T. Yakhshiev, Eur. Phys. J. A36 (2008) 37 

45. J. Gasser and H. Leutwyler, Phys. Rept. 87 (1982) 77 

46. R.G Thomas, Phys. Rev. 81(1951)148 

47. J.B. Ehrman, Phys. Rev. 81 (1951) 412 

48. A. Bulgac, V.R. Shaginyan, Nucl. Phys. A 601 (1996) 
103 

49. S.A. Fayans, JETPLett. 68 (1998) 169 

50. A. Bulgac, V.R. Shaginyan, Eur. Phys. J. A 5 (1999) 
247 

51. S.A. Fayans and D. Zawischa, Int. J. Mod. Phys. B 15 
(2001) 1684 

52. B.A. Brown, Phys. Rev. C 58 (1998) 220 

53. B.A. Brown, RIKEN Rev. 26 (2000) 53 



